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Abstract. A general phenomenological reaction-diffusion model for flow-induced phase transitions in com- 
plex fluids is presented. The model consists of an equation of motion for a nonconserved composition vari- 
able, coupled to a Newtonian stress relations for the reactant and product species. Multivalued reaction 
terms allow for different homogeneous phases to coexist with each other, resulting in banded composition 
and shear rate profiles. The one-dimensional equation of motion is evolved from a random initial state to 
its final steady-state. We find that the system chooses banded states over homogeneous states, depending 
on the shape of the stress constitutive curve and the magnitude of the diffusion coefficient. Banding in the 
flow gradient direction under shear rate control is observed for shear-thinning transitions, while banding 
in the vorticity direction under stress control is observed for shear-thickening transitions. 

PACS. 47.20. Ft Instability of shear flows - 47.20.Hw Fluid dynamics: Morphological instability; phase 
changes - 05.45.-a Nonlinear dynamics and nonlinear dynamic systems - 05.70.Ln Nonequilibrium and 
irreversible thermodynamics 



1 Introduction 

There is a significant body of experimental evidence doc- 
umenting the existence of sharp, stable interfaces sepa- 
rating two or more phases or "bands" , in shear flow in 
complex fluids. This phenomena has been reported in var- 
ious types of surfactant solutions jj} ^polymers ||, liquid 
crystals j| and colloidal suspensions ||. There appears to 
be a compelling generality between these "phase transi- 
tions" in different complex fluids: 

i) The onset of banding or phase separation manifests 
itself as a discontinuity in the "flow curve" of the sys- 
tem. The flow curve is the unique relationship between 
the measured shear stress and the applied shear rate 
(or vice versa) at steady-state. An experimental flow 
curve typically contains segments that correspond to 
homogeneous flow, as well as segments corresponding 
to inhomogeneous flow. The individual homogeneous 
bands which make up the inhomogeneous state each 
have their own homogeneous flow curve, which we shall 
refer to as a "constitutive curve" . The inhomogeneous 
flow curve then represents the response of the system, 
averaged over different spatial regions that occupy dif- 
ferent homogeneous flow branches, in proportions to 
maintain the externally controlled shear stress or shear 
rate. [Henceforth we will use the terms stress and shear 
stress interchangeably, unless otherwise specified.] 



ii) The transition only occurs above a unique and repro- 
ducible critical stress or shear rate. 

iii) The flow curve can be qualitatively different depend- 
ing on whether the average stress or the average shear 
rate in the system is held fixed. [In a typical rheo- 
logical experiment this is achieved by controlling the 
torque or angular velocity respectively] For intermedi- 
ate stresses or shear rates, the flow curve usually has 
multiple branches which are not equally accessible un- 
der both stress and shear rate control. For weak and 
strong flows, the flow curve is single-valued, and the 
same locus of points is traced out under stress or shear 
rate control. 

iv) The flow-induced bands have different shear rates or 
shear stresses, and are generally also distinguished by 
some combination of different degrees of order and dif- 
ferent microstructures. 

v) The interfaces between the bands may be aligned in 
the direction of the flow gradient or the flow vortic- 
ity. Each banding orientation has its own rheological 
signature. In shear-thinning systems, for example, a 
stress plateau in the flow curve usually indicates gradi- 
ent banding, while extrema in the stress (as a function 
of shear rate) usually indicate vorticity banding. One 
of us H has constructed possible flow curves based 
on the banding orientation and the character of phase 
coexistence (shear-thinning versus shear-thickening). 
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Gradient banding has been unambiguously observed in 
solutions of wormlike micelles. In strain-controlled exper- 
iments on shear-thinning solutions, a stress plateau coin- 
cides with shear-banding in the gradient direction H. In 
shear-thickening solutions H], a gel-like phase can be in- 
duced by flow. Under shear rate control the induced phase 
fills the system at steady-state, resulting in a discontinu- 
ous stress jump in the flow curve. Under stress control 
phase coexistence between solution and gel is observed, 
the gel fraction being an increasing function of stress. In 
the corresponding flow curve the shear rate shows a min- 
imum and maximum. 

Vorticity banding has been reported in dense colloidal 
suspensions || and surfactant solutions of multilamellar 
vesicles When the shear rate is held fixed, the flow 
curve shows a maximum and minimum in the stress. Un- 
der controlled stress, there is a jump up in shear rate 
upon increasing stress, and a jump down in shear rate 
upon decreasing stress. The same qualitative curves have 
also been observed in surfactant hexagonal phases || , al- 
though in that case vorticity banding has not yet been ex- 
plicitly verified. Such behavior is analogous to that of the 
shear-thickening wormlike micelles, if the roles of stress 
and shear rate are interchanged. In bcc cubic crystals 
of triblock copolymers Eiser et al. observe two stress 
plateaus in the flow curve under controlled shear rate. X- 
ray diffraction shows that each plateau corresponds to dif- 
ferent orientations (relative to the flow direction) of dense 
planes in the crystal. 

In steady-state there can be no acceleration, so the to- 
tal stress must be divergence free. In planar shear flow, this 
implies that the shear rate in the vorticity direction and 
the shear stress in the gradient direction are uniform. Vor- 
ticity banding thus corresponds to a scenario where bands 
share a common shear rate but can have different shear 
stresses (see Figure |l|) . Similarly, when bands lie in the gra- 
dient direction the stress is uniform across the bands and 
the shear rate can vary. Most experiments where banding 
has been observed have been carried out in the curved 
geometries of cone-and-plate or Couette rheometers. The 
gaps in these rheometers are usually very thin, and in this 
limit the flow is approximately a planar shear flow. [We 
also note that we consider flows in the low Reynolds num- 
ber limit.] 

The microscopic mechanisms causing the transitions 
in all these complex fluids are likely to be highly system 
specific, and govern the critical shear rate or shear stress, 
the structure of the flow-induced phases, and the detailed 
shape of the flow curves. At a macroscopic level, however, 
there appears to be a high degree of universality between 
systems. As we have discussed, different complex fluids 
can produce qualitatively similar flow curves. By simply 
analysing the shape of these flow curves, we have extracted 
information about the banding orientation of the transi- 
tion H , as well as the stability of the system || . 

The obvious analog to this way of thinking is the well- 
known Landau- Ginzburg theory of equilibrium phase tran- 
sitions. A free energy functional consisting of a double- well 
local free energy and a square gradient term reproduces 




Fig. 1. a) Gradient banding: flow- induced phases lie in the 
direction of the velocity gradient (arrow shows flow direction). 
The bands share the same shear stress, but can have different 
shear rates, b) Vorticity banding: bands lie in the flow vorticity 
direction. Here, the shear rate is uniform across the bands, but 
the shear stress can vary from one band to the other. 

all the phenomenology of a phase transition in the region 
of the critical temperature, for many different systems. 
However, a microscopic theory is required to calculate the 
Landau coefficients. In this work, we use a multivalued re- 
action diffusion scheme to construct a general phenomeno- 
logical theory to describe phase transitions in flow; in the 
spirit of Landau- Ginzburg theory, such a model could in 
principle be derived from microscopic theories. 

The most difficult step in devising such a nonequilib- 
rium theory is determining the important variables. Typ- 
ically there are three significant quantities: i) a species 
concentration, which is a conserved scalar, ii) the mo- 
mentum density, which is a conserved vector and iii) the 
stress, which is a non-conserved tensor. The momentum 
density is described by the Navier-Stokes equation, and 
its current is the stress. In addition non-conserved "mi- 
crostructural" order parameters exist which contribute to 
the stress. These may be scalars such as chain length in 
wormlike micelle solutions, or tensors such as molecule ori- 
entation. Since all of these variables have different relative 
relaxation times, one must distinguish between slow vari- 
ables, which require their own equations of motion, and 
fast variables, which relax quickly to a steady-state value. 

The choice of slow variables affects the structure of the 
equations of motion and the couplings between them, and 
therefore the dynamics of the system. In models of hydro- 
dynamic instabilities for example, the momentum is con- 
sidered to be a slow variable. A phenomenological stress 
constitutive equation is often used: if the stress is taken 
to be a fast variable, this relation is simply an algebraic 
function of the rate of strain tensor, such as the Newtonian 
relation for simple fluids; if the stress is taken to be a slow 
variable, this relation takes the form of a differential equa- 
tion, such as the Upper Convected/Oldroyd-B Maxwell 
model for polymer melts. The hallmark of complex fluid 
rheology, however, is the coupling between the velocity 
and/or the stress to the microstructure of the fluid. In 
microscopic theories, generally an equation of motion is 
not written for the total stress, but for another slow vari- 
able which makes an important contribution to it, such 
as the director in nematic liquid crystals or the second 
moment of the configuration tensor in polymer melts. 
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Schmitt, Marques and Lequeux |lCj have classified flow 
instabilities in complex fluids as "mechanical" or "spin- 
odal" instabilities, using a model where concentration and 
momentum are the slow variables. If a perturbation to the 
shear rate first makes the system go unstable, the instabil- 
ity is mechanical, while it is spinodal if the concentration 
becomes unstable first. Note that "instability" as it is used 
here refers to a linear instability. Any instability, linear or 
nonlinear (we return to this issue at the end of the paper), 
can lead to a macroscopically shear-banded state that re- 
solves the instability. 

Shear-banding associated with momentum instabilities 
have been analyzed in detail at a high (macroscopic) level, 
using the phenomenological Johnson-Segalman model Ji"l[ 
|l^,|l3|. Here a non-conserved "polymer" stress tensor, play- 
ing the role of the slow variable, is included in the total 
momentum density, resulting in a multi-valued stress con- 
stitutive relation. This model produces gradient banding 
and a flow curve with a stress plateau, and is consid- 
ered a reasonable mimic of shear-thinning wormlike mi- 
celles. Microscopically derived theories for wormlikc mi- 
celles [[l4| and nematic liquid crystalline melts |25| yield a 
non-monotonic relation similar to the Johnson-Segalman 
model, but with the benefit of a molecular interpretation. 

An alternative caricature to these models has been de- 
veloped in phenomenological theories for shear-thickening. 
Originally, Ajdari Jl5| proposed an equation of motion for 
the position of an interface that separates high and low 
viscosity phases under shear. By coupling this equation 
with conservation laws and a Newtonian stress constitu- 
tive equation for the micellar solution, a non-monotonic 
flow curve was produced. Goveas and Pine adopted this 
approach to describe shear-thickening wormlike micelles 
and were able to successfully reproduce much of the ex- 
perimental phenomenology. The flow curve was then used 
to explain the differences in stress versus shear rate con- 
trol, based on a linear stability analysis of the interfacial 
height equation. In this case, the momentum and micel- 
lar solution stress were taken to be fast variables, while 
an equation of motion was written for a scalar variable, 
which is the macroscopic manifestation of changes in the 
fluid microstructure. 

However, the formulation of Goveas and Pine did not 
contain any mechanism for the formation of the shear- 
induced state, so that the existence of the new phase was 
simply postulated by the presence of an interface. In this 
paper we present a generic phenomenological model that 
naturally admits a flow-induced phase and incorporates 
spatial gradients so that an interface structure and its sta- 
bility can be determined. The model consists of an equa- 
tion of motion for the volume fraction of a reacting species, 
i.e. a scalar non-conserved order parameter representing 
microstructural change in a complex fluid. There are fast 
stress variables associated with the reactant and product 
species, which contribute additivcly to the total stress. We 
have continuously evolved the model from a homogeneous 
to a phase-separated state, and examined how thinning 
and thickening flow curves, as well as the size of the gra- 
dient terms, affect phase transitions, and in particular the 



banding orientation (vorticity versus gradient banding). 
Most significantly, we are able to probe the nonlinear dy- 
namic behavior of the system. 



2 Minimal model 

Our phenomenological theory consists of a general reaction- 
diffusion scheme. The reaction terms represent the cre- 
ation and destruction of a variable under flow, and are 
analogous to the local free energy terms in a Landau- 
Ginzburg theory. This variable may embody a species con- 
centration, or a structural parameter such as aggregate 
size or molecule orientation. While this scheme is meant 
to be quite general and is a vehicle for capturing the 
general physics for many complex fluids, a reaction dif- 
fusion scheme has a literal basis for wormlike micelles and 
onion solutions. In the wormlike micelle case, such "re- 
action" terms correspond to the constant breaking and 
recombination of the "living" polymers; while in onion 
solutions, the reaction terms might correspond to the for- 
mation of onions. The steady-state onion size scales as the 
inverse square root of the shear rate jL(| and is a reversible 
function of the shear rate; i.e. the size is independent of 
whether smaller onions are created by increasing the shear 
rate applied to larger onions, or larger onions are created 
by decreasing the shear rate applied to smaller onions. 
This indicates that onion combination and fracture pro- 
cesses compete to attain steady state, and these processes 
have different dependences on shear rate. 

In the same way that a double well potential signals 
the possibility of equilibrium phase coexistence, a multi- 
valued reaction term can allow for flow-induced banding. 
The diffusion terms are the analog of the non-local terms 
in the free energy and provide gradients which can support 
inhomogeneities and describe interfaces between states. 
However, unlike in equilibrium, where a global minimiza- 
tion principle applies, the diffusion terms are necessary for 
determining the conditions for phase coexistence in flow 

s 

Consider a system which is one-phase at equilibrium, 
and consists solely of a species, A. Planar shear flow is then 
applied to this system: the coordinate system is shown in 
Figure |], where x, y and z denote the flow, gradient and 
vorticity directions respectively. We consider only varia- 
tions in y and z in this work. Suppose that a new phase, 
B, can be induced by flow, such that at a given shear 
rate (or stress) a dynamic equilibrium between A and B 
is established. We write this schematically as: 

A^B. (2.1) 

Let us define an order parameter, <pB = 4>i which corre- 
sponds to the volume fraction of the B-species (although 
allowing cj) to correspond to a structural variable is also vi- 
able). The system is constrained to have constant density 
such that 

(j> A + <p B = l- (2.2) 
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Notice that 4>a and <f>B are non-conserved variables, al- 
though the total density is conserved. We write the fol- 
lowing equation of motion for as 



(2.3) 



where i?(0, 7) represents the forward and backward "reac- 
tions" which create and destroy the new phase, j(y, z) = 
dv x /dy is the local shear rate (where v x is the component 
of the velocity in the flow direction), and D is an effec- 
tive diffusion coefficient (taken to be a consta nt). The ho- 
mogeneous steady-state solutions to equation ( |2.3| ), where 
V 2 — 0, are given by 



i?(^,7)=0. 



(2.4) 



3 Constitutive curves 



To compute flow curves for the system, stress constitu- 
tive equations for the different components must also be 
specified. The simplest possible scheme involves additive 
Newtonian relations for each species 



(7 = OA + (7B 

0~a = -qa^aj, 



(3.1) 



where a is the total shear stress, and o~ a and r\ a are the 
shear stress and viscosity, respectively, of species a = 
A, B. Stress and composition are thus effectively coupled 
i n th e system. Applying the density constraint, equation 
(2.2), gives an expression for the constitutive curve of the 
system, 



<7=[0(C-1) + 1]7 > 



(3.2) 



where = 0(7) is the solution to equation (2^), c = 
rjB /tja is the ratio of the viscosities of the two components, 
and we have set t\a = 1 for simplicity. 

Notice that a multivalued i?(0, 7) produces multiple 
homogeneous steady-states. Choosing 



(3.3) 



yields the curve shown in Figure g, for an imposed 7. [k 
represents a rate constant for a backward reaction, which 
has dimensions of inverse time and is henceforth set to 
unity .] N otice that the forward reaction term in equa- 
tion ( |3.3| ) has a linear shear rate dependence, so we must 
take its absolute value from symmetry considerations. In 
microscopic theories for flow-induced reactions in worm- 
like micelles |l7j] and polymers Jl8| , such an effective reac- 
tion rate also has a linear or non-analytic form, resulting 
from the projection of tensorial degrees of freedom onto a 
scalar order parameter. For a given local shear rate, the 
reaction scheme of equations (2.2, 2.4, |3.3|) yields the fol- 



lowing homogeneous steady-states: 
0!=O 



1 1 L 4 

63 =2 + 2V^7- 



(3.4a) 
(3.4b) 

(3.4c) 
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Fig. 2. Multivalued reaction scheme (see equation 2.3) pro- 
ducing three possible roots, 4>i,4>^ an d 03, corresponding to 
stable, unstable and stable homogeneous states, respectively, 
for a fixed value of the shear rate. 



Performing a linear stability analysis on equation (2.3), at 
fixed 7, shows that 4>i and </>3 are stable fixed points for 
the system, while <p2 is an unstable fixed point. This is 
evident from Figure g, where the middle root has positive 
slope dR/dcj) > 0. 

Substituting the homogeneous roots from equation ([TJ) 
into equation ( |3.2| ) produces three branches of the consti- 
tutive curve (labelled as 1,2 and 3 respectively in Figure 
||). Notice that below a certain shear rate, 7 = 4, only the 
4>\ root is real and the reaction curve is single- valued. This 
is marked as point P on Figure ||, and has coordinates 



{7p,^p} = {4,2(c+1)}. 



(3.5) 



Physically, this means that only specie s A exists at low 
(uniform) shear rates. From equation (3.2), we can see 
that this one-phase system (4>a = 1) is Newtonian, and the 
corresponding flow curve has a slope of unity. The slope 
of the stable flow-induced branch 3 depends on the value 
of the parameter c. For c < 1, a transition from branch 
1 to branch 3 is shear-thinning, while for c > 1 such a 
transition is shear-thickening. The locus of flow induced 
roots (02,03) exhibits a minimum in the shear stress as 
a function of shear rate, which is denoted as point Q in 
Figure 0: 



{jq,o-q} = I 



(1 + ^) 5 



,(1 + V~c)' 



(3.6) 
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Fig. 3. Flow curves (from equations 2.4, 
sponding to the minimal model. At low shear rates, the flow 
curve is single- valued (branch f ). Above a certain shear rate (or 
stress), an additional unstable branch, 2, and a stable branch, 
3, exist. These branches correspond to the homogeneous roots 
{(f>i,4>2,4>i\ from equations (1.4) for imposed shear rate, and 
{(f>i, 4>2, 03 } from equations (3.7) for imposed stress. In the for- 
mer case, branches 2 and 3 are separated by P, while they are 
separated by Q in the latter. Thus, the line segment between 
P and Q is stable under controlled shear rate, but unstable un- 
der controlled stress, (a) Shear-thickening flow curve for c > 1, 
illustrating the controlled shear rate case. At a fixed shear 
rate, P, the system can choose between homogeneous states 
on branches 1 and 3, or gradient band between these branches 
at stress a*, (b) Shear-thinning flow curve for c < 1, illustrat- 
ing the constrolled stress case. At fixed stress E, the system 
can vorticity band at shear rate 7* between 1 and 3, or choose 
between homogeneous states on these same branches. 



For a given value of the local stress, the homogeneous 
steady-states are given by 



k = 

(o--c+l) 



»2' = 



4cr 



2(7 



(3.7a) 
(3.7b) 



2cr 



Note that these roots are found by recasting equation fl3.3|) 
in terms of the stress, by using equation (3^2). This proce- 



dure is no t eq uivalent to substituting equation (3.2) into 
equation ([0|). This is because while the loci of homoge- 
neous states is the same under fixed local stress or shear 
rate, the stability of these steady-states is not; i.e. the por- 
tion of the constitutive curve between P and Q is unstable 
under fixed local stress, but stable under fixed local shear 
rate. In Figure ||, point Q marks the stress above which 
the constitutive curve is multivalued for controlled stress, 
while P marks the strain rate above which the constitutive 
curve is multivalued for controlled shear rate. 

In an experiment, however, only the average stress and 
shear rate can be controlled. If the average shear rate is 
held fixed at P, for example, the system can choose be- 
tween various options (illustrated in Figure ^a): 

i) A homogeneous low stress state, (pi. 

ii) A homogeneous high stress state, 4>y . 

iii) A mixture of states (i) and (ii), where the interfaces 
between phases lies in the vorticity direction (vorticity 
banding). Note that (i) and (ii) cannot coexist with 
each other in the y-direction, since the stress must be 
homogeneous in the gradient direction. 

iv) A mixture of high shear rate phase, <j>y & n d a low shear 
rate phase, <pi. Here the system attains an intermedi- 
ate stress, (J*, and the relative proportions of the two 
phases are set such that the average shear rate is main- 
tained at r. Since the bands have the same stress, but 
different shear rates, this scenario corresponds to gra- 
dient banding. 

While it appears from Figure ^ that there is a multiplicity 
of stresses a* at which the system can gradient band, in 
fact the system selects a particular stress (see next Sec- 
tion) . 

If instead, the average stress is fixed at S (illustrated 
in Figure |^b) , the system can choose a high or low shear 
rate homogeneous phase, or it can gradient band between 
these. Alternatively, it can band in the vorticity direction 
between high and low stress states, by adopting a shear 
rate 7*. There is also a selected shear rate for vorticity 
banding. The essential question is: which of the many pos- 
sible states available to it does the system actually choose 
and why? 



4 Calculating the banding stress and shear 
rate 



The inclusion of gradient terms in equation (2.3) causes 
stress selection for gradient banding, and shear rate se- 
lection for vorticity banding jt9). The selected stress and 
shear rate are determined by mathematically connecting 
two different homogeneous stable states to form an inho- 
mogeneous profile. To find the banding shea r ra te 7* at 
which vorticity banding can occur, equation ( |2.3| ) is inte- 
grated across the domain at steady-state. A banding solu- 
tion (homogeneous phases separated by interfaces) is, by 
definition, one which has no gradients in ip at the bound- 
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aries. We obtain the following condition: 

r<fe(7*) 



d(j)R[<j)^*] = 0. 



(4.1) 



01 (7*) 



where cf>i and ^3 are given by equations ( |3.4| ). Defining 
a new fun ctio n, _F[(/>,7] = d<j)'R[(j)' , 7], we can rewrite 
equation ( fP] ) as 



(4.2) 



which is analogous to the common tangent construction 
from equilibrium thermodynamics. If this were an equi- 
librium system, F would be identified as the free energy. 
Since this is a dynamic system however, F cannot be given 
the same physical interpretati on, so that the anal ogy is 
purely formal. Using equation fl3.4|) in equation (4.2) gives 



7* = 4.5. 



(4.3) 



To calculate the banding stre ss a* at which gradient 
banding can occur, equation (4.1) must be recast in terms 
of the she ar s tress using the stress constitutive relation, 
equation (3.2) to obtain a relation 7(0", cj>). The banding 
stress in our minimal model is only a function of c and is 
given by the solution of the following equation: 



F[4> v (a*),a*]=F[4> 3 ,(a*),a*], 
where F[<p,a*} = d(j)' R[(j>' ,7(0-, <p')] yields 



(4.4) 



{c-iy 

- (2c +1)4> 



3 + i(c+2)0 2 + |c (4.5) 



113 
3 



2 

cln 



(4.6) 



and 4>v and <p3' are given by equations ( |3.7| ) . The selected 
stress a* and shear rate 7* are shown in Figure ^ for c = 
0.3, 2.7, and given in Table [| 



c 


* 




7(^3) 


op 


7Q 


a Q 


0.3 


2.815 


0.810 


6.504 


2.6 


4.374 


2.395 


0.6 


3.606 


0.732 


5.100 


3.2 


4.066 


3.149 


1.2 


4.910 


0.642 


4.352 


2.4 


4.008 


4.391 


2.7 


7.625 


0.529 


4.014 


5.4 


4.252 


6.986 



Table 1. Banding stress a*, points of instability P and Q, 
and coexistence conditions for different values of c. In all cases 
7P = 4 and 7* = 4.5., while the stress and shear rate on branch 
(f>i are related by <j\ = 71 . 



5 Dynamical Selection of Steady-States 

In the preceding Sections we have seen that certain ho- 
mogeneous and banded states are available to the system, 



based on a steady-state analysis. To determin e wh ich of 
these states is selected in practice, equation (2.2) must 
be evolved in time; to make contact with experiments we 
can only impose constraints of fixed average shear rate or 
stress. In this work only one-dimensional calculations are 
performed, so that the equation of motion is solved either 
in the y (gradient) or z (vorticity) directions. If there are 
composition modulations in the y-direction, these can only 
cause modulations in the shear rate (gradient banding), 
since the shear stress must be uniform in y. If the average 
stress is controlled, there is only one "interesting" stress 
a* at which the system can become inhomogeneous in the 
y direction. Due to the numerical difficulty of fixing a pre- 
cise stress in the system, we do not consider this case. By 
comparison, if the average shear rate is controlled, there 
is a wide range of shear rates for which we can investi- 
gate whether the system remains homogeneous or gradi- 
ent bands. Similarly, when composition modulations in z 
are allowed, we can only look for vorticity banding under 
controlled stress within this calculation. Gradient banding 
can occur only if the shear rate is set exactly at 7*, which 
we do not study here. 



5.1 Controlled Average Shear Rate 

We first consider the system under shear rate control, 
and only allow for spatial var iations in y. Integrating the 
stress relation, equation (3.2), across the domain (where 
the shear stress is independent of y) gives the local shear 
rate as a function of V, the velocity difference across the 
system. Then equation ( |2.3| ) becomes 



V 



^ 2 (W) 



dt f dy {!/[</> (c-l) + l}} 



D 



dy 2 ' 
(5.1) 



This is an integro-differential equation, instead of the dif- 
ferential equation which yielded the analysis of Section |^. 
Notice, however, that the sam e homogeneous steady-states 
are obtained. Equation (5.1) is solved using random ini- 
tial conditions with <f> uniformly chosen within the range 
[0 — 1], and no flux boundary conditions, keeping V at 
a fixed value. The domain size is normalized to unity, so 
that V is synonymous with the average shear rate. 

To solve equations (5.1), we use a fully implicit finite 
difference scheme, using a central difference approxima- 
tion for first and second spatial derivatives, and a forward 
difference approximation for the time derivative. Nonlin- 
ear terms are linearized in time as follows: 



W[(j>{x, t + At)} = W[<j>{t)} + [<j>{t + At) - <j>(t)} 



dW[4>,t] 

d4> ' 

(5.2) 



The integral in equation (5.1) is evaluated explictly, i.e. at 
the previous time step. In general, 300 spatial mesh points 
are used with a time-step of 1/10000 (2yj. 

For some initial conditions, the resulting steady-states 
are homogeneous, while for others, they are banded. This 
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Fig. 4. Basins of attraction for different states, for c = 0.3 
and imposed mean shear rates V = 7. In I the system usually 
attains the homogeneous flow branch = <f>i = 0. In II the 
system usually bands at stress a* , and in III the system usually 
attains the homogeneous state <f> = 03 . The behavior is smooth 
as V is increased through these regimes (see data in Table ^). 

implies that there is a basin of attraction for attaining 
a banded state. Our results can be categorized accord- 
ing to the shape of the constitutive curve and the magni- 
tude of the diffusion coefficient, and whether the system 
is shear-thinning or thickening. We find that decreasing 
the diffusion coefficient increases the basin of attraction of 
the banded state and destabilizes the homogeneous state. 
That is, the system is more likely to band for narrow inter- 
faces. Intuitively, this makes sense since the banded state 
represents a mathematical connection between two homo- 
geneous states: the wider the interface, the more difficult 
it is for gradients to be non-zero near the boundaries of the 
system. For c < 1, for V < 7p the system usually chooses 
branch 1 of the constitutive curve over the banded state, 
even for small diffusion coefficients. Notice that banding is 
first allowed, in principle, when the imposed shear rate is 
larger than that of the low shear rate band. The "critical" 
shear rate, where the system actually first starts to band, 
is generally somewhat higher than this: it is the low shear 
rate limit of the stress plateau that would be measured 
experimentally. Below the critical shear rate the system 
always chooses branch 1, although the exact location of 
this point depends on the diffusion coefficient. [Increas- 
ing the diffusion coefficient widens the interface, affecting 
where the interface first "touches" the wall [|l3| and the 
ability of a banded system to satisfy the boundary condi- 
tions, as discussed above.] 

When 7p < V < 7q, the system is more inclined to 
band rather than remain homogeneous, even for large dif- 
fusion coefficients (see Fig. [|). The chosen homogeneous 
states always lie on branch 3 of the constitutive curve. Be- 
yond 7q, there is a crossover to where homogeneous states 
(on branch 3) are preferred over banded states. The exact 
crossover depends on the value of the diffusion coefficient. 
This shows that the system has made a transition under 
flow, so that for c < I a shear-thinning transition is seen, 
as discussed in Section |§|. These findings are illustrated in 



Table 2, for 49 runs with different initial conditions with 
c = 0.3. The "critical" shear rate in this case is around 
V = 3.3. Figures H and H show similar behavior for c = 0.6. 

Gradient banding is never observed in our numerical 
experiments when c > 1. Here, the chosen final steady 
state, above the critical shear rate, is always the stable 
high stress homogeneous state on branch 3, which makes 
this a shear-thickening transition. 



Table 2. Summary of results for c = 0.3, for 49 runs with 
random initial conditions and controlled mean shear rate V. 
In this case, jp corresponds to V = 4.0, and 7q corresponds 
to V = 4.374. Banding is first allowed at V = 2.815. 
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The preceding results are for the case where the shear 
rate is held at a steady value, and might apply to a sys- 
tem where the mean shear rate is applied to an initially 
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noisy system. However, most experiments are conducted 
by starting up the system from zero shear rate, and then 
discontinuously ramping the shear rate to higher values. 
We have tried to mimic this scenario, by bringing the sys- 
tem to steady-state for V < jp, and then suddenly in- 
creasing V to greater than jp. In order to dislodge the 
system from branch 1 to the banded state, however, we 
need to add noise of amplitude order unity. Such a large 
amount of noise essentially obliterates any memory of the 
initial steady-state, suggesting that a nucleation event is 
required for an experimental system to band from start- 
up, as was found previously in Ref. [ fL3| . 

In general, the banded state consists of two bands, cor- 
responding to branches 1 and 3 on the constitutive curve. 
For some initial conditions, multiple bands are found. We 
do not attach much significance to this result, because it 
is known that multiple interfaces are allowed for planar 
flow, and the number of allowed interfaces in such one- 
dimensional systems is known to increase as the diffusion 
coefficient decreases Britton and Callaghan have re- 
ported multiple gradient bands for wormlike micelles in 
Couette flow fjj. However, it has been shown jl3| that 
simple constitutive relations (like the one derived here) 
do not permit multiple interfaces in Couette flow. This 
implies that the current model cannot describe these ex- 
perimental observations. 



5.2 Controlled Average Stress 

Next we fix the average shear stress and solve the model 
in the z-direction, hence allowing for vorticity banding at 
different stresses. Following a procedure analogous to that 
of the previous Section, the observation that t he s hear rate 
is uniform in z allows us to convert equation (2.3) into the 
following integro-differential equation 



27</> 2 (l - 



dt fdz[(/>(c-l) + l] 



n-£, (5.3) 



where 27 is the imposed average shear stress across the 
domain. 

Note that interchanging the stress and shear rate vari- 
ables in the phenomenological model would "reverse" all 
the results of Section 5.1 . Thus, branch 3 would be thicken- 



ing for c < 1 , and the system would be inclined to vorticity 
band under stress control, for op < 27 < <jq. Similarly, 
the system would remain homogeneous for c < 1. This 
reasoning implies that the shape of the flow curve sets the 
attractors for the system. If this is true, then when stress 
and shear rate are not interchanged, the system should be 
predisposed to vorticity band for c > 1 when the stress is 
fixed between ap and ctq, but should remain homogeneous 
for all stresses when c < 1. 

This is ind eed what we find from numerical solution of 
equation ( [3.3] ) . Table || shows results for various values of 
27 and c = 2.7. Decreasing the diffusion coefficient desta- 
bilizes the homogeneous state, as in the controlled shear 
rate case. 



V=4.05, C=0.6, D=0.01 



Volume Fraction as a function of Spatial Position 



V=4.05, C=0.6, D=0.005 








Volume Fraction as a function of Spatial Position 
Fig. 5. Steady-state composition profiles for 49 runs with 

dom initial conditions, for c = 0.6, V = 4.05. Here, 7p cc 

sponds to V = 4.0, and 7q = 4.066. 
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V=4.1,C=0.6, D=0.001 










Volume Fraction as a function of Spatial Position 
V=4.1, C=0.6, D=0.01 



Volume Fraction as a function of Spatial Position 
Fig. 6. Steady-state composition profiles for 49 runs with ran- 
dom initial conditions, for c = 0.6, V = 4.1 



Table 3. Summary of results for c = 2.7, for 49 runs with 
random initial conditions and controlled stress S. In this case, 
ap = 5.4 and <jq = 6.896. 
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We sometimes observe multiple bands under stress con- 
trol, as with the fixed shear rate cases. Bonn et al. 0] 
and Chen et al. have seen multiple bands in the vortic- 
ity direction in Couette flow which is probably due to 
a combination of the slow coarsening expected in one- 
dimensional systems | p2| and multiple allowed interfaces 
The stress is non-uniform and monotonic in the flow 
gradient direction of a cylindrical Couette device, which 
implies a single stable interface. The cylindrical geometry 
does not, however, impose such an inhomogeneity along 
the vorticity direction. 



6 Conclusions 

We have shown that a simple phenomenological reaction- 
diffusion scheme can produce a flow-induced phase tran- 
sition, as a consequence of a multi-valued reaction term. 
The model consists of an equation of motion for a non- 
conserved composition variable, while the stresses induced 
in the reactants and products are assumed to be fast vari- 
ables. The character of the model depends on a single pa- 
rameter c, that controls whether or not the transition is 
shear-thinning or shear-thickening. Above a critical shear 
rate (or shear stress), the system may band or remain 
homogeneous. The steady-states that are selected from 
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random initial conditions depend on the shape of the con- 
stitutive curves and the magnitude of the diffusion coeffi- 
cient: 

1. Imposed shear rates: For c < 1 (shear-thinning tran- 
sition), the system chooses a low stress homogeneous 
state at low shear rates. Above a critical shear rate, 
gradient banding tends to occur for imposed shear 
rates around the region of the constitutive curve with 
negative slope da/dj < (see Figure ||a). At shear 
rates higher than this, the system is predisposed to- 
wards the high stress homogeneous state. For c > 1 
(shear-thickening transition), the system always chooses 
this homogeneous state above the critical shear rate 
and gradient banding is never observed. 

2. Imposed stress: For c > 1 (shear-thickening transition) , 
the system chooses a low shear rate homogeneous at 
low stresses. Above a critical stress, vorticity band- 
ing tends to occur for imposed stresses around the re- 
gion of the constitutive curve with negative slope. For 
higher stresses, the system is predisposed towards the 
high shear rate homogeneous state. For c < 1, the sys- 
tem always chooses this homogeneous state above the 
critical stress and vorticity banding is never observed. 

3. In the regions of the flow curve where banding is ob- 
served, we find the apparent basin of attraction for 
banding increases upon decreasing the value of the dif- 
fusion coefficient. 

While banding is more pronounced in the vicinity of the 
flow curve with a negative slope (where the system is lin- 
early unstable), it is also observed in regions of the flow 
curve with positive slope. In particular, the critical shear 
rate or stress (where banding is first initiated), lies in the 
latter section of the flow curve. Here, the system is nonlin- 
early unstable to perturbations. Such behavior has been 
seen in experiments on shear-thinning wormlikc micelles 
p3f , where the onset of banding occurs at a lower stress 
(and shear rate) if the system is given enough time to 
explore all fluctuations, as compared to where banding 
is induced upon rapidly varying the control parameters. 
Porte et al. [24j have discussed various flow curves which 
can contain both linearly and nonlinearly unstable regions: 
the equilibrium analog of the former is the spinodal curve, 
and that of the latter is the metastable region, where an 
instability must be nucleated. 

Our results are significant because they show that a 
minimal model can exhibit a rich phenomenology, and 
that the selection rules for phase coexistence are simple. 
To understand why the system chooses certain states over 
others in some regions of the flow curve, a nonlinear dy- 
namics analysis of the model must be performed. We be- 
lieve that our scheme represents a new class of reaction- 
diffusion equations, because the constraint of fixed average 
stress or shear rate turns the governing partial differential 
equations into integro-differential equations, which repre- 
sents a general class of dynamical equations that, to our 
knowledge, has not been studied. This system exhibits fas- 
cinating and complex nonlinear dynamics, which we will 
discuss in a future publication. 



Our current scheme is missing much physics: a com- 
plete model would involve coupled equations of motion for 
conserved variables (concentration of the various species) 
and non-conserved, tensorial variables (structural variables, 
stress). Also, we have assumed that the individual species 
obey Newtonian stress constitutive relations. Typically, 
these species are themselves complex fluids, and are either 
shear-thinning or exhibit a yield stress. In future work, 
by systematic exclusion of certain dynamic variables, we 
will be able to investigate the individual roles played by 
the stress, concentration etc., in order to determine which 
variables are essential to the problem formulation. 

One of us p5[ has already considered a theory with 
stress, momentum and concentration variables in the con- 
text of rigid rod suspensions. Separate phase diagrams for 
shear-induced phase separation in both the vorticity and 
gradient directions were calculated, but the model was 
too prohibitively complicated to study which of these ori- 
entations would in fact be selected by the system. In this 
work, we have used a much simpler scheme to demonstrate 
the neccessary analysis (albeit within a one-dimensional 
model- see next paragraph) to unambiguously determine 
whether banding actually occurs in a system, as well as 
the banding orientation. While Schmitt et al. have 
also presented quite a simple phenomenological model (in- 
cluding both concentration and momentum as dynamical 
variables), they did not go beyond a linear stability anal- 
ysis. They also did not consider the case (as we have here) 
of a non-conserved variable initiating an instability in the 
system. 

Our calculations have been carried out only for the 
case of planar flow. It has been shown for the Johnson- 
Segalman model |l3| that the nonuniformity of stress in 
a curved geometry has significant effects on banding. In 
addition, we have examined the issue of gradient versus 
vorticity banding using a one-dimensional model. Realis- 
tically, the model should be solved considering both vor- 
ticity and gradient directions simultaneously. The band 
orientation may be influenced by anisotropy in the dif- 
fusion coefficient. A convective term of the form v • V^> 
should also be included in the equation of motion. Such a 
term does not appear in a one-dimensional shear flow, but 
it can qualitatively affect transients in phase separation in 
two dimensions. Finally, noise has been incorporated into 
our model through the initial conditions. While Gaussian 
noise is present in the equation of motion through the dif- 
fusion term, in a driven system there may be other noise 
terms that should be added. 

We stress that our phenomenological theory only aims 
to describe the general macroscopic physics of flow-induced 
phase transitions. Details concerning the underlying struc- 
tural transformations can only be probed by more specific 
microscopic models. 
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